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Abstract 



O . Relativistic Mean Field calculations have been performed for a number 

' of nuclei in mass A ~80 region. Ground state binding energy, charge 

radius and charge density values have been compared with experiment. 
Optical potential have been generated folding the nuclear density with the 
microscopic nuclear interaction DDM3Y. S-factors for low energy (p, 7) 
and (p,n) reactions have been calculated and compared with experiment. 

Relativistic Mean Field (RMF) approach has proved to be very successful 
in explaining different features of stable and exotic nuclei like ground state 
^ . binding energy, deformation, radius, excited states, spin-orbit splitting, neutron 

f-H ' halo, etc|T]. Particularly, the radius and the nuclear density are known to be 

well reproduced. Its has led to its application to nuclear reactions also. 

Low energy reactions are very important from the astrophysical point of 
view. In astrophysical environments, neutron and proton reactions are the keys 
QQ ' to nucleosynthesis of heavy elements. Mass 80 region is an interesting one as 

00 . {p, 7) and {n,p) reactions play important roles in determining the abundance of 

' elements. In this mass region, there are some proton-rich naturally occurring 

isotopes of elements known as p nuclei which can not be produced via s- or 
r- process. Mainly, proton capture reactions contribute to the formation of 
such nuclei. Recent works [See e.g. [2]] have emphasized the importance of a 
number of charge exchange reactions in this mass region for production of very 
light p nuclei. The relevant astrophysical rates can directly be derived from the 
{p, n) data, though the target is in the ground state and reaction has negative 
, a Q-value[3]. A number of recent experiments has focussed on reactions using 

^ ' protons having a few MeV energy in mass 80 region. 

Calculation of isotopic abundance requires a network calculation involving 
many thousands of reactions. Despite the importance of (p, 7) or {p, n) reactions 
in explaining the abundance of p nuclei, experimental data are rather scarce 
due unavailability of the target nuclei on earth. Thus, one often has to depend 
on theory for these reactions. The calculations essentially utilize the Hauser- 
Feshbach formalism where, the optical model potential, a key ingredient, is 
often taken in a local or a global form. It is also possible to use a microscopic 
optical potential constructed utilizing nuclear densities. If the target is stable, 
the density of the nucleus is available through electron scattering. However, in 
absence of stable target, theory remains our sole guide to describe the density. 
Thus, it is imperative to test the theoretical calculations, where experimental 
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data are available, to verify its applicability. We aim to to check the success 
of microscopic optical potentials based on mean field densities in explaining 
the reaction cross sections. A good description depending essentially on theory 
will allow one to extend the present method to the critical reactions which are 
beyond present day laboratory capabilities. 

Calculations using microscopic potentials have been able to explain the ob- 
served elastic scattering cross sections even in nuclei far off the stability valley 
(See, e.g. Ref. H] and references therein). Low energy projectiles probe only 
the outermost regions of the target nuclei. Hence, the nuclear skin plays a very 
important role in such reactions. The density information should be available 
from theoretical calculations. This method has been utilized to study low en- 
ergy proton capture reactions in Ni and Cu nuclei [5_ and nuclei in A 60 — 80 
region [6]. 

For the present study we have selected a number of low energy proton re- 
actions for their astrophysical relevance. The reactions ®'*'*^'®^Sr(p, 7) were in- 
vestigated through activation technique by Gyiirky et al. in Ref. [7]. It is 
important to note that *'*Sr is another p nucleus. In beam measurements were 
performed by Galanopoulos et aL[5] to find out the cross sections for the re- 
action ^^Sr(j), 7). As for charge exchange reactions, three reactions, identified 
as important by Rapp et a/. [5] and for which experimental cross sections are 
available, have been selected for study. The reaction '^^As(p, n) was studied 
in [21 [TUJ [TT] through in beam detection of neutrons. Finally, the reactions 
'^^Ge(p, n)[T2] and *^Rb(p, n)[3] were studied by activation technique. In the 
present work, we investigate the reactions mentioned above in a microscopic 
approach. 

Theoretical density profiles were extracted in the RMF approach. There are 
different variations of the Lagrangian density as well as a number of different 
parametrizations. In the present work we have employed the FSU Gold|13j 
Lagrangian density. It contains, apart from the usual terms for a nucleon me- 
son system, nonlinear terms involving self coupling of scalar-isoscalar meson, 
and additional terms describing self-coupling of the vector-isoscalar meson and 
coupling between the vector-isoscalar meson and the vector-isovector meson. 

Pairing has been introduced under the BCS approximation using a zero 
range pairing force of strength 300 MeV-/m for both proton and neutrons. 
The RMF-(-BCS equations are solved under the usual assumptions of classical 
meson fields, time reversal symmetry, no-sea contribution, etc. Since we need 
the densities in co-ordinate space, the Dirac and the Klein Gordon equations 
have directly been solved in that space. This approach has earlier been used 
[H im [TS] in neutron rich nuclei in different mass regions. 

The microscopic optical model potentials for the reactions are obtained us- 
ing effective interactions derived from the nuclear matter calculation in the lo- 
cal density approximation, i.e. by substituting the nuclear matter density with 
the density distribution of the finite nucleus. In the present work, the micro- 
scopic nuclear potentials have been constructed by folding the density dependent 
DDM3Y|16l [T7] effective interaction with the densities from the RMF calcula- 
tion. This interaction, obtained from a finite range energy independent M3Y 
interaction by adding a zero range energy dependent pseudopotential and intro- 
ducing a density dependent factor, has been employed successfully in nucleon 
nucleus as well as nucleus nucleus scattering, calculation of proton radioactivity, 
etc. The density dependence has been chosen in the form C(l — (3p^^^)\XI\, the 
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constants being obtained from a nuclear matter calculation [TB]. The real and 
the imaginary parts of the potential are taken as 0.7 times and 0.1 times the 
DDM3Y potential, respectively. This normalization have also been used in our 
earlier work on (^,7) reactions in lighter nuclei[B]. We have checked that the 
above values adequately describe the cross section measurements. Of course, 
these parameters can be tuned to fit the cross sections in individual reactions. 
For example, in ^^Sr, if we choose the imaginary part of the potential as 0.3 
times of the DDM3Y potential, the result will differ by 10% and fit the experi- 
mental data better. However, we believe a single parametrization for the entire 
mass region to be more useful. 

The Coulomb potentials are similarly constructed by folding the Coulomb 
interaction with the microscopic proton densities. We have already used such 
potentials to calculate life times for proton, alpha and cluster radioactivitv[T^ 
as well as elastic proton scattering^ in different mass regions of the periodic 
table. 

Reaction calculations have been performed with the computer code TALYS 
1.2[5^ assuming spherical symmetry for the target nuclei. DDM3Y interaction 
is not a standard input of TALYS but can easily be incorporated. Though the 
nuclear matter-nucleon potential does not include a spin-orbit term, the code 
provides a spin-orbit potential from the Scheerbaum prescription [2 II coupled 
with the phenomenological complex potential depths. The default form for this 
potential given in the code has been used without any modification. 

The TALYS code has a number of other useful features. We have employed 
the full Hauser-Feshbach calculation with transmission coefficients averaged over 
total angular momentum values and with corrections due to width fluctuations. 
Up to thirty discrete levels of the nuclei involved have been included in the 
calculation. 

Our calculations, being more microscopic, are more restricting. Yet, the rate 
depends on the models of the level density and the El gamma strength function 
adopted in the calculation of cross sections. Phenomenological models are usu- 
ally fine tuned for nuclei near the stability valley. Microscopic prescriptions, on 
the other hand, can be extended to the drip lines, and hence, have been assumed 
in all nuclei. We have calculated our results with microscopic level densities in 
Hartree-Fock (HF) and Hartree-Fock-Bogoliubov (HFB) methods, calculated in 
TALYS by Goriley and Hilaire, respectively. We have also compared our re- 
sults using phenomenological level densities from constant temperature Fermi 
gas model, back shifted Fermi gas model and generalised superfluid model from 
TALYS. The cross sections are very much dependent on the level density chosen, 
sometimes changing by a factor of 50%. We find that in most of the cases, the 
HFB densities fit the experimental data better in our formalism. 

For El gamma strength functions, results derived from HF-I-BCS and HFB 
calculations, available in the TALYS data base, have been employed. In agree- 
ment with our observation in in [6], here also the results for HFB calculations 
describe the S-factors reasonably well and we present our results for that ap- 
proach only. 

It is possible to scale the theoretical capture cross sections to match with 
experiment using a parameter Gnorm in the code used to scale the gamma-ray 
transmission coefficient. However, for the present paper, we have not scaled 
the theoretical results. All the parameters in the Lagrangian density and the 
interaction are standard ones and have not been changed. 
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As the density profile is the important quantity in our formalism, a com- 
parison of the radii values can provide some idea about the agreement of the 
calculated densities with experiments. In Table 1, we compare our results for 
the binding energy and charge radii (rch) with measurements for those nuclei in 
this mass region, which have been used as targets for low energy proton capture 
or charge exchange reactions. The binding energy values from the mean field 
approach have been corrected using the formalism developed in [S^ . The 
experimental binding energy values are from Ref. |24| . 

Charge radii have been calculated from the charge densities, which, in turn, 
have been obtained from the calculated point proton density pp by taking into 
account the finite size of the proton. The point proton density is convoluted 
with a Gaussian form factor g(r), 

Pchir)^ J epp{v')g{v-v')dv' (1) 
5(r) = (a0F)-3exp(-rVa') (2) 

with a — 0.8 fm. 

Experimental charge radii values are from Angeli[25]. The results show that 
RMF can describe the charge radii of these nuclei with sufficient accuracy. One 
sees that in most of the nuclei, the difference between measurement and theory 
is less than 1%. 

Direct comparison of charge density is more difficult in absence of accurate 
experimental information. De Vries et a/. [55] have presented the coefficients 
of Fourier-Bessel expansion for charge density of a number on nuclei extracted 
from electron scattering data. It includes two nuclei of our interest, ''^Ge and 
^*Sr. In Fig. [U we compare the charge density extracted from the Fourier- 
Bessel coefficients and our calculated results for the above two nuclei. One can 
see that the theoretical and experimental values agree very well, particularly 
at larger radii values, which is the region expected to contribute to the optical 
potential at low projectile energy. However, in absence of information on error 
in the density values, this conclusion can remain only tentative. 

Next, we compare the results for the reaction calculation in the above men- 
tioned reactions with experiments. As the astrophysically important Gamow 
window lies in the region 1.3 to 3.9 MeV for these nuclei, we present the re- 
sults covering this energy region. The cross-section varies very rapidly at such 
low energy making comparison between theory and experiment rather difficult. 
The usual practice in low energy nuclear reaction is to compare another key ob- 
servable, viz. S-factor. The expression of the astrophysical S-factor^ is given 

by, 

S{E) = Ecj{E)e^''''' (3) 

where E is the energy in centre of mass frame in KeV, cr{E) is reaction cross- 
section in barn and 77 indicates the Sommerfeld parameter which may be ob- 
tained from the relation, 27rr; = 'il.2QZpZt\J p./ E. Here, Zp and Zt are the 
charge numbers of the projectile and the target, respectively and p is the re- 
duced mass (in amu). It varies much slowly than reaction cross-sections as the 
exponential energy dependence of cross-section is not present in it. For this rea- 
son, we calculate this quantity and compare it with experimentally extracted 
values. 
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Figures 2 and 3 show the results for the reactions ®^'^^~®*Sr(p, 7). The 
resuhs compare favourably with experiments compared to the NON-SMOKER 
code calculations of Rauscher et al.^7\. However, it needs to be pointed out 
that, in the case of ®''Sr, theoretical results overpredict the cross section values. 
It was suggested fj] that perhaps the agreement with theory (in their case the 
NON-SMOKER calculation) worsens as one goes to more neutron-rich nuclei. 
However, as one can see in the right panel of Figure 2, this trend is not shared 
by the present calculation. 

Figure 3 shows the results for the {p,n) reactions on (a) ^^As, (b) ^^Ge and 
(c) ^^Rb targets. These reactions (along with their inverse reactions) are listed 
among the ten most important reactions in deciding the abundance of the p 
nuclei in Rapp et al.^. The three measurements for the ''^As{p,n) reaction 
are rather old and error values are not available for most of the measurements. 
The quoted error in cross section is 10% or above. In ^^Ge, we find that the 
calculation systematically overpredict the results by as much as 60%. On the 
other hand, the calculations for the *^Rb(p, n) reaction produce excellent match 
with experimental measurements. 

We find that our calculation can reproduce the S-factor values with reason- 
able success. Even in the worst case, calculation is off by a factor less than two 
while the cross section values range over four orders of magnitude. However, one 
should remember that in astrophysical calculations, the rates are often varied by 
a large factor, viz. ten or hundred[5]. Thus the present microscopic calculations 
can be used to obtain rates which are dependable for astrophysical calculations. 

We point out that in our earlier work [5 , we have seen that the default local 
and global optical potentials [5S] in the TALYS package also can be used with 
suitable normalization of gamma ray strength to produce comparable results 
for certain energy ranges. In the present case also, suitable selection of the 
parameter brings the values calculated with default potential close to experi- 
mental values. However, we believe that the present microscopic approach is 
more suitable as no normalization is necessary and the method can be extended 
to reactions where experimental data are not available. 

In summary, cross sections for low energy (p, 7) and {p, n) reactions for a 
number of nuclei in mass 80 region in the energy regime important for explo- 
sive nucleosynthesis have been calculated using the TALYS code. The micro- 
scopic optical potential has been obtained by folding the DDM3Y microscopic 
interaction with the nuclear densities obtained from RMF calculation using the 
Lagrangian density FSU Gold. 
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sored DRS Programme of the Department of Physics of the University of Cal- 
cutta. CL acknowledges the grant of a fellowship awarded by the UGC. GG 
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Table 1: Experimental binding energy and charge radii values compared with 
calculated results. 





B.E.(MeV) 


Teh 


(fm) 




Exp. 


Theo. 


Exp. 


Theo. 




728.90 


727.53 


4.236 


4.232 




748.93 


748.27 


4.226 


4.240 


«7Sr 


757.36 


757.17 


4.220 


4.245 


88sr 


768.47 


768.47* 


4.220 


4.249 


75 As 


652.56 


652.38 


4.097 


4.082 


76 Ge 


661.60 


660.69 


4.081 


4.053 


85Rb 


739.28 


738.70 


4.203 


4.218 



* Normalized following the prescription of [22l [23] 
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Figure 1: Comparison of charge density obtained from Fourier-Bessel analysis of 
experimental electron scattering data (solid line) and calculated in the present 
work (dashed line). 
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Figure 4: Experimental and calculated S-factors for (a) ^^As{p,n), (b) 
'^^Ge{p,n) and (c) *^Rb(p, n) reactions, respectively. 
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